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Abstract 

A generalization of the Kohn-Sham approach is derived where the correlation-energy functional 
depends on the one-particle density matrix of noninteracting states and on the external potential 
from the interacting target-state. The one-particle equations contain the exact exchange potential, 
a nonlocal correlation potential, and an additional operator involving the correlation density. The 
electronic-energy functional has multiple solutions: Any one-particle density matrix delivering the 
target-state density yields a solution. In order to obtain the Kohn-Sham solution, the nonlocal 
operators are converted into local ones using an approach developed by Sala and Gorling. Since the 
exact exchange-potential is used, and the A^-represent ability problem does not arise — in contrast 
to the Kohn-Sham approach — errors from Coulomb self- interactions do not occur, nor the need to 
introduce functionals defined by a constraint search. Furthermore, the approach does not use the 
Hohenberg-Kohn theorem. A density functional formalism is also derived that assumes that the 
one-particle density matrices of interest have v-representable (non-interacting) densities and that 
these density matrices can be written as an explicit functional of the electron density. For simplicity, 
we only consider noninteracting closed-shell states and target states that are nondegenerate, singlet 
ground-states. 



I. INTRODUCTION 



The Kohn-Sham version of density functional theor 



chemistry and condensed matter physics 




[lays a major role in both quantum 
,|7|. The local density approximation 
has been widely used for the solid state. While for molecules, by far, the most successful 
functional, a hybrid one ^, [l^ U) 1^; is known as B3LYP jol [l^. 

The Kohn-Sham approach, however, does have well known shortcomings. For example, a 



constraint search definition 
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17( 1 is required to treat the f-representability problem 



that arises in the original Kohn-Sham method j^. Unfortunately, this formal definition is 
difficult to consider when deriving approximate functionals. Furthermore, in contrast to 
wave function based methods, the exchange- correlation functional is an unknown, implicit 
functional, and there is no systematic method to improve approximations. In addition, 
there are well known errors ari sing from Coulomb self-interactions that appears when using 
approximate functionals 0, Q|. Also, the most widely used approximate functional for 
molecular systems, the B3LYP functional, includes a component of the exact exchange- 

;he noninteracting state to come 
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25| is an 



potential, even though the Kohn-Sham approach requires 
from a local potential. The optimized potential method 
approach to convert a nonlocal operator into a local potential. Unfortunately, this method 
leads to potentials that are not invariant to a unitary transformation of orbitals and depend 
explicitly on the individual orbitals and orbital energies. 

The formalism presented below uses an electronic-energy functional containing a correla- 
tion energy functional Sco that depends on the external potential v and on the one-particle 
density matrix pi of determinantal states. Since the f-representability problem does not 
appear, a constrain search definition is not needed. Also, since the approach uses the ex- 
act exchange-potential, errors from Coulomb self-interactions do not occur. The energy 
functionals, however, contains multiple solutions, since any one-particle density matrix pi 
delivering the density from the interacting state yields a solution. In order to obtain the 
Kohn-Sham solution, the nonlocal operators are converted into local ones using an ap- 
roach developed by Sala and Gorling j^]. In contrast to the optimized potential method 



Hacn aeveiopea D 
0,0, HQ, 



25|, the energy functionals and local potentials are invariant to a 



unitary transformation of orbitals and do not depend on the individual orbital or the orbital 
energies. A density functional formalism is also derived that assumes that the one-particle 
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density matrices of interest have ti-representable (non-interacting) densities and that these 
density matrices can be written as an exphcit functional of the electron density. 
Previous 



theory 



y we have shown that the correlation energy from many b ody perturbation 
3, 28,0 can be written as an explicit functional of v and pi j30|. In a similar 



manner, but using less restrictive energy denominators, the correlation energy functionals 
presented below can be shown to be an explicit functional of v and pi 0|. Hence, in contrast 
to the Kohn-Sham method, it maybe possible to derive approximate functionals that can 
be improved in a systematic manner. For simplicity, we only consider noninteracting closed- 
shell states and target states that are nondegenerate, singlet ground-states. 

II. THE ENERGY FUNCTIONALS AND TRIAL WAVE FUNCTIONS 

Our interest is in finding the ground-state eigenvalue of the Hamiltonian operator, 

H^, = f + V,, + V,, (1) 

where 

N 

T = E(-i^')' (2) 



1 ^ 

V.. = ^rT-, (3) 

N 

K = (4) 

i 

and V is the external potential; is the number of electrons. Since the Hamiltonian Hnv 
is determined by and f, so are the ground state wave functions \^ nv) that satisfy the 
Schrodinger equation: 

Hnv\^Nv) = £nv\^Nv)-, (5) 

where, for simplicity, we only consider wave functions that are nondegenerate, singlet ground- 
states. 

Using a second quantization approach, our spin-free Hamiltonian does not depend on N , 
and it can be expressed by 

H. = 5^(z|(-iV2)|j)4- + J2^Mj)E,, + 1 J2itj\kl)E,,,i, (6) 

ij ij ijkl 



where the symmetry-adapted excitation operators are given by 

Eijki = '^al^al^aacija, (8) 

crA 

and the one- and two electrons integrals are spin-free integrals written in chemist's notation 
js^ using a spatial orbital set, say {x}] this set has the following form: 

^ia(x) = Xi(r)cr(a;); cT = a,/3, (9) 

where the spatial and spin coordinates, r and u, are denoted collectively by x. 

Wave function-based methods including perturbation theory, configuration interaction, 
and coupled cluster theory, use one or more reference states to express \E' and S. For closed- 
shell ground-state wave functions, a single determinant can be used, where closed-shell 
determinant al, or noninteracting, states can be constructed from a set of doubly occupied 
spatial-orbitals; these occupied orbitals also determine the spin-less one-particle density- 
matrix of the noninteracting state, given by 



Pi(ri,r2) = 2 J2 X-(ri)x:(r2), (10) 

where the sum is over the occupied orbitals; this set of orbitals is denoted by {Xo}- 
For later use, we also mention that for a complete basis set we have 

2(5(ri - r2) = pi(ri,r2) + Kp,(ri,r2), (11) 

where is determined by the excited orbitals, 

Kp,(ri,r2) = 2 J2 Xr(ri)x:(r2), (12) 

and {Xu} denotes the set of orbitals orthogonal to the occupied set {Xo}- The operator form 
of Eq. (HID is 

2i = pi + kp„ (13) 

where I is the identity operator; so, the kernels of the three operators within Eq. p3p are 
given by the corresponding terms within Eq. (jllll . 



It is well known that there is a one-to-one mapping between determinantal states and 



their one-particle density matrices [2, l34l, say 7, where for a closed-shell state described by 
the orbitals given by Eq. 0, we have 35|, |36|, |37|, |38[ 

7(xi,X2)= J^x^(ri)x^(r2)CT(a;i)a*(^2), (14) 

"'e{xo} o" 

and by using Eq. (jTU)), we obtain 

7(xi,X2) = ^Pi(ri,r2)(5<^i<^2. (15) 
Since our closed-shell determinantal states are determined by pi, we denote these kets by |pi). 



According to the Hohenberg-Kohn theorem ^^Q;!^' external potential v is deter- 
mined by the density, and the density also determines A^. So, in principle, we can replace the 
variables N and v by the electronic density n and, at least for nondegenerate ground-states, 
write 

H,\^n)=£n\'^n); n-^ N,V, (16) 

where these functions serve as density-dependent trial-wave functions for the Kohn-Sham 
approach. Notice we have omitted the N subscript on the Hamiltonian operator, since Hy 
is independent of when this operator is expressed in second quantization. 

As an alternative to a density- dependent wave function, we consider trial wave functions, 
say |\&t,pi), that are determined by the one-body external potential v and, in addition, by 
the spin-less one-particle density-matrix pi of a noninteracting state, and, as mentioned 
previously, these noninteracting states are denoted by \pi). 

By definition, our trial wave function {"^vpi) yields the exact ground-state wave function 
when the noninteracting density ps, i.e., the density of |pi), equals the exact density n 
of the interacting state l^&n), where n also determines the v and A^. This state of affairs can 
be represented by the following: 

\^vp,) = l^n); Pi — > Ps = n, n — ^ A, V. (17) 

In other words, pi determines ps, and when ps = n, \'^vpi) yields l^^n)- Letting gi denote 
the one-particle density matrix of interest, we can write 

\^v^,) = \^n); Qi — >ri, n — >N,v. (18) 
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For later use, we also mention that the density n of an interacting state can be partitioned 

as 



n = ps + Pc, (19) 



where the correlation density is given by 



, , (^„|f(r)|^„) 



and r is the density operator, given by Eq. ()A7|) . 

Using our trial wave function, we introduce a variational energy functional: 

EM - 1^ . • (21) 

Our trial wave functions \^vpi) and energy functionals are assumed to be explicit 

functionals of pi and v. However, two different one-particle density matrices, say pi and 
Pi', that yield the same density p^, i.e., pi — > ps and pi — > ps, yield the same \'^vpi) 
and E^[pi], so these functions are implicit functionals of Ps, and, therefore, we can write 
l^^vps) and Ey[ps]. However, we will continue to consider them as functionals of their explicit 
variable pi. 

Using Eqs. ffT^ and (fTHjl. we observe that our energy functional E^, given by Eq. (j^H), 
delivers the exact energy Sn when the one-particle density matrix determines the exact 
density n: 

EM]=^n, Qi — ^n, n — ^A^,^, (22) 
and for an arbitrary density we get 

E,[pi]>£n, pi-^Ps^N, (23) 
where the density ps from the noninteracting state \pi) is not necessarily f-representable. 

III. TRIAL HAMILTONIANS 

Our trial wave function is a ground-state eigenfunction of a Hamiltonian operator that 
depend explicitly on the one-particle density of a noninteracting state: 

H^pA^vp,) = SvpA^vp^)- (24) 



As in our trial wave functions I^I/tjpi) and energy functionals E^[pi\, the trial Hamiltonians 
Hvp^ are explicit functionals of pi, but implicit functionals of ps- So two trial Hamiltonians, 
say Hvp^ and Hvpi', are equal if both pi and pi yield the same density, i.e., pi, pi ps- 

Our trial Hamiltonians must be chosen so that Eq. (fTH|l is satisfied, indicating the follow- 
ing identity: 

H^g^ = Hy, Qi — > n, n — ^A^,^;. (25) 

There are many ways to obtain a trial Hamiltonian that satisfies Eq. fj25p . Consider the 
following trial Hamiltonian obtained by adding a term to the Hamiltonian: 

H^p^ = + \ J drgpXr) (^f (r) - p,(r)j , pi — > p^, (26) 
where r(r) is the density operator, given by Eq. ()A7|1 : (f (r) — Ps(rV) is the one-body portion 



H 0, Q 



given by 



of r(r) when this operator is written in normal-ordered form 
Eq. ()A6|) . Furthermore, A is an arbitrary constant, and the functional g is also arbitrary, 
except that it vanishes when the correlation density pc vanishes 

lim«,^(r) = 0, (27) 

where pc is defined by Eqs. (fT9|) and (f20|) . 

Since (r(r) — Ps(r)) is normal-ordered, we have 

(pi|(f(r)-p,(r))|pi)=0. (28) 

Therefore, the added term appearing in Eq. ()2(i|l can be considered a sort of correlation 
term, since it does not contribute in first order. Hence, we have 

{pi\H,p,\pi) = {pi\H,\pi). (29) 

One possible choice for Qp^, and presented in Appendix |X1 is given by 

gpX^i) = j dY2r^2 Pc{tc2)- (30) 

IV. A GENERALIZATION OF THE KOHN-SHAM FORMALISM 

We now obtain a generalization of the Kohn-Sham formalism. Substituting Eq. ^ into 
Eq. dn} gives 







^.Pl) 







E.[p,] = + ^"^^^-'^^^X^ + / dvv{v)p,{v) + I rfr.(r)p.(r), (31) 



where pc is the correlation density of the trial wave function, i.e, as in Eq. fl2Up . we have 

p,(r) = i^%^Htf^-p,(r)=n-p.(r), ^ n, Pi ^ p., (32) 

and n is the density of "^vp^- 

Through the first-order, the kinetic energy and electron-electron repulsion energy are 
given, respectively, by 

(pilTlpi) = j dv, [-|V?Pi(ri,r2)]r^^^^, (33) 
(pilKelPi) = Ej[p,]+E^[p^], (34) 
where the Coulomb and exchange energies are 

Ej[ps] = \jj ri2dridY2p{Yi) p{y2), (35) 
1 f f 



-E^pi] = - J J rfrirfr2pi(ri,r2)pi(r2,ri). (36) 
Adding and subtracting (pi|T|pi) and {pi\Vee\pi), Eq. (jHU can be written as 
EM = j dr, [-iV?pi(ri,r2)]^^^^^ + J rfrt;(r)p,(r) 

+ Ej[ps] + E4pi] + E^4pi,v] + J drvir)p^{r), (37) 
where the correlation-energy functional is given by 

Recognizing the first four terms from Eq. (jHTj) as the energy through the first order. Si, we 
can write 

E^[pi] = £i[pi,v]+ Eco[pi,v] + j drv{r)pc{r), (39) 

where 

Si[pi,v] = {pi\HM = J dri [-|V?Pi(ri>r2)]r,=r, (40) 

+ J driv{ri)p{ri) + ^ J J dridv2r^2 p{ri) p{r2) - ^ j j dridr2r\2Pii^i^'^2)piir2,ri). 

Now consider the correlation energy that is obtained by wave function methods. Using 
the notation from Eq. (jSJ, and a reference state |pi), the correlation energy is given by 

Sco[pi,v]= Si[pi,v], (41) 



where previously we have shown that £co can be written as an exphcit functional of v and 
pi js^. In a similar manner, but using less restrictive energy denominators, our correlation 
energy functional Eco-, given by Eq. (j38p . can be shown to be an explicit functional of v and 
Pi |3l|. Therefore, by requiring the last term within Eq. (j39j) to be an exp licit functional of 
V and pi, can also be written as an explicit functional of v and pi j31[. 

We now focus our attention on minimizing the energy functional subject to the 
constraint that the spin-less one-particle density-matrix pi comes from a closed-shell single- 
determinantal state. For the more general case of a determinantal state, say I7), with the 
(spin-dependent) one-particle density matrix 7, as in Eq. (fT^ . the two necessary conditions 
for 7 to satisfy are given by the following 3, Q : 

7(X3, X4)5(X3 - X4) (iX3dX4 = A^, (42) 
j 7(X3, X5)7(X5, X4) (ix5 = 7(X3, X4), (43) 

where the first relation indicates that the electron density yields the number of electrons A^; 
the second relation indicates that 7 is indempotent. For our special closed-shell case, we 
substitute Eq. (fT3j) into the above constrains, yielding the following conditions: 




Pi(r3,r4)(5(r3 - r4) <ir3<ir4 = A^, (44) 
j pi(r3, r5)pi(r5, r4) drg = 2pi(r3, r4). (45) 

It is well know that the functional derivative of £1 with respect to the 7 yields the kernel 
of the Fock operator JJ]. For the closed-shell case, we have 

i^ri,r2 =:r^ X) (46) 

'5pi(r2,ri) 

where, using Eq. (001), the Fock kernel is given by 

Fp,(ri,r2) = 5(ri - V2) + v{r2) + j dv^r:,^ p{v-i)^ +<Hri,r2), (47) 

and the exchange operator, say v^} , has the following kernel: 

<nri,r2) = -^rr2Vi(ri,r2). (48) 

By generalizing Eq. ()4fi|l . we define a generalized, or exact, Fock operator JF, where the 
kernel of this operator is 

^,,(ri,r2) = = F,,(ri,r2) +<Hri,r2) +C(ri,r2), (49) 



and the correlation operator v^^ and external-correlation operator -O^^^ are defined by their 
kernels: 

<c 1-1,^2 = . . ^ • (51) 

Minimizing the functional E^, given by Eq. subject to the constraints given by 

Eqs. (j44j) and (|45p . is very similar to the corresponding Hartree-Fock derivation |^ and 
the derivation for reference-state one-particle density matrix theory The only 



difference being that the spin variable has been eliminated, and we have a factor of two 
appearing in Eq. ()45|) . Therefore, we only state the main results, i.e., this minimization 
yields the exact electronic energy Sn for the interacting state, as given by Eq. (j22j) . where 
the one-particle density-matrix satisfies the following conditions: 

ke^Pq^Qi = 0, (52) 

and the kernels of the operators pi and /tp^ are given by the terms on the right side of 
Eq. dH}; also, as mentioned previously, gi yields the exact density n of the interacting 
state Using Eqs. (jH^ and (jSSI), it is readily shown that jF^j and commute: 



0, (54) 



and the occupied orbitals satisfy a generalized Hartree-Fock equation: 

oiXui ^ ^ ^xw~X.x-i (55) 

where the notation x & qi indicates a summation over the occupied orbitals from the deter- 
minantal state Xw is also an occupied orbital from Furthermore, we can choose 
orbitals that diagonalize the matrix e^,^, yielding exact, canonical Hartree-Fock equations: 

{-\V^ + v + v] + «f + vil + v',l ) Xw = e^Xw , Xw^Qu (56) 

where the Coulomb operator is defined by 

^^j(ri)x(ri) = j rfr2rf2V(r2)x(ri), (57) 
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and we have 

^i(r,r) =n(r). (58) 

Equation ()56|) is also satisfied by the canonical excited orbitals. 

For later use, we also mention that the determinantal states \gi) satisfy the following 
noninteracting Schrodinger equation: 

5^^,,(r,)|^i)=2(5^eJ |^>i). (59) 

Appendix IbI presents an alternative way of partitioning the energy functional that differs 
from Eq. 

V. CONVERSION OF THE NONLOCAL POTENTIAL INTO A LOCAL ONE 

As mentioned previously, our energy functionals are implicit functionals of the non- 
interacting density Ps- Hence, any one-particle density- matrix that yields the interacting 
density minimizes our energy functional, i.e., we have 

Sn = [pi] = E, \q\] = E, [g'l] ■ ■ ■ , (60) 

where 

n(r) = ^i(r, r) = ^^(r, r) = ^"(r, r) ■ ■ ■ , (61) 

and there are other solutions besides Eq. (jHHjl . e.g, 

^q[Xw = (-|V^ + v + vj + w^'J Xw = e'ujXw, Xw e Qi, (62) 

where the nonlocal potential Wp^ is given by 

Wp,=K'+v',l+v^l. (63) 

Assuming n is a noninteracting v-representable density, there exist a noninteracting state, 
say that has n as its density: 

n(r) = V5i(r,r), (64) 
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and this determinant — assuming it is a closed-shell determinant — is the ground-state solution 
of the following noninteracting Schrodinger equation: 

f;/(r,)|(/.i) = 2(5]ej \^,), (65) 

i=l \ w J 



where 



f = -lV' + Vs, (66) 



and Vs is a local potential. Therefore, the canonical occupied orbitals from satisfy the 
following one-particle Schrodinger equation: 

f(l)w={-\^^ + V + vJ + V^c)(pw = (^w4>w, e (67) 

where with no loss of generality, we have required Vg to be defined by 

Vs = v + v^ + v^c- (68) 

By definition, or using Eqs. (jHIIl), (j^ - and ^pi is a one-particle density matrix that 
minimizes our energy functional: 

Sn = E4^i\, (69) 
and, therefore, cpi also satisfies Eq. (j^ : 

f;^^,(rO|^i) = 2[^eJ l^i). (70) 



i=l 



Hence, it follows from Eqs. ()65p and ()7Up that \(fi) is an eigenstate of two different nonin- 
teracting Hamiltonians. By comparing Eq. (jU^ and (jU7j) with g[ = ipi, we see that the two 
operators, JF^^ and /, are identical, except that J^^-^ contains the nonlocal operator w^-^ and 
/ contains the local potential v^c- Furthermore, the occupied orbitals from Eq. ()fi2j) and ()fi7|l 
with q[ = Lpi may differ by a unitary transformation, but they yield the same one-particle 
density matrix: 

M^i,r2) = 2 X-(ri)x;(r2) = 2 ^ ^^1)^(^2). (71) 



Using the approach by Sala and Gorling [26|, and Eqs. (jHSI), (fTOj) . (jHH) and (jU7j) . but 
permitting the orbitals to be complex, it is readily demonstrated that v^c is given by 

^^xnfr) = y dri [2«;(ri,r)(/?i(r,ri) - (/?i(r, ri) J rfrg v?i(r2, r)ty(ri, rs) (72) 



-h(/?i(ri,r)v3i(r,ri)fxc(ri) 
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By substituting v^c repeatedly on the right side we can obtain an expansion for v^c- 
[2«;(ri, r)v9i(r, ri) - V9i(r, ri)(y9i(r2, r)w(ri, r2) 



2n(r 

+ V5i(ri,r)v9i(r,ri 



—^{w{r2, ri)v?i(ri, rz) - ^¥?i(ri, r2)(/Ji(r3, ri)«;(r2, rg)} 
n(ri) 2 



+ (/7i(ri,r)(/?i(r,r 



1 



-V?i(r2,ri)v9i(ri,r2) 



w(r3,r2)(/7i(r2,r3) + 



(73) 



VV 1 L X llXl uXlV^ 

j; is the Kohn-Sham 



where there are integrations over the dummy variables ri, r2 and r^. The leading term of 
Eq. dZSl) is the Slater potential HHQ; 

this term also appears within the Krieger-Li- 
lafrate (KLI) approximation of the optimized potential method 
The orbitals (pw satisfying Eq. ()fj7|l are the Kohn-Sham orbitals 
noninteracting state. However, / differs from the Kohn-Sham operator, since, in addition to 
depending explicitly ipi, instead of n, / depends explicitly on the external potential v from 
the interacting Hamiltonian H^j. Furthermore, the external-correlation operator vj^^ does not 
appear in Kohn-Sham formalism. In addition, unlike the original Kohn-Sham approach , 
the iV-representability problem does not arise, nor the need to introduce a constraint-search 



definition 
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3,0, 



17l | to avoid this problem. 



In our derivation we have assumed that \lpi) is a ground state solution of Eq. ()65p. 
However, the results may also be valid if \ipi) is an excited state solution, since the Sala and 
Gorling approach may also be valid in this case. 



VI. CONVERSION OF THE ONE-PARTICLE DENSITY-MATRIX FUNCTION- 
ALS INTO DENSITY FUNCTIONALS 

For noninteracting states, the wave function is determined by the one-particle density ma- 
trix. For certain closed-shell determinantal states, we can write pi [ps] , where this functional 
includes all densities that are noninteracting f-representable, but it is also defined for all 

nnnn 

A^-representable densities. Using the constraint search approach [IJ, ll5|, [la III! , for a given 



density, say p', the functional pi[p'] yields the one-particle density matrix that minimizes 
the expectation value of the kinetic energy: 



Min,(pi|T|pi) = (pi[p']|T|pi[p']), 
pi ^ p 



(74) 



where the search is over all determinantal states that have a density of p'. 
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Substituting pi [p] into Eco of Eq. (|H7|) gives 

E,[p,] = j dn [-lV?pi(ri,r2)]r^^r^ + j dTv{v)ps{r) 

+ Ej[p,]+E,[p] + E,,[p,v] + j drv{r)p,{r), p ^ pi, (75) 

where, using pi[p], the last term is also a functional of v and p. This equation differs 
from the Kohn-Shan density functional, since the correlation-energy functional depends on 
the external potential v, and the last term does not appear in the Kohn-Sham approach. 
However, mathematically speaking, the minimization of Eq. ()75|) follows the same procedure 
as in the Kohn-Sham method, yielding 

Uw= {-^"^^ + v + v^ + v^ + v^, + v^,)(t)^ = e^(l)^, 4>we^i, (76) 

where the local potentials are given by 

^P(r) = '^^-[^'^J (77) 
^xirj ^^^^^ , {II) 

. . ^ 5E^^[p,v] 
^coirj ^^^^^ , {1^) 

^ecirj - ^^^^^ . 

Assuming the density n from the interacting state is noninteracting z;-representable, we 
have 

Ey[n] = £n, n is noninteracting i;-representable. (80) 

Note that Eq. (f75j) is a valid energy functional only when the one-particle density matrix 
that enters the first term is the same one generated by the functional pi[p]; this is the case, 
at least when p is non-interacting t^-representable. 

APPENDIX A: A POSSIBLE CHOICE FOR Qp^ 

The electron-electron repulsion operator is spin-free and can be written as 

Vee = \Y.mr^2\km3ku (Al) 
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where the two-electron integral is written in chemist's notation [32| and the two-electron 
spin-adapted excitation-operator is given by Eq. (jH)). This operator can also be written as 



where the pair-function operator is given by 

f2(r2,ri) = i^Xi(ri)x*(ri)xKr2)Xfc(r2)^yfci, 



(A2) 



(A3) 



ijkl 



and this operator yields the diagonal elements of the spinless two-particle density matrix as 



the expectation value. Writing this operator in normal-ordered form |27l . 
respect to the vacuum state |pi), we have 

= 



40 



41 



m with 






dridr2r-^^2 P2{r2,ri)p^ + J J dridr2r-^2 Psir2)T{ri 

dridr2r^2 Pii^2,ri)t{ri,r2)p^ + I / dridr2ri2^f2(r2,ri)p^,(A4) 



(A5) 



where, examining each term in turn, from the first term we have 

1 



P2(r2,ri)p, = -p,(r2)p,(ri) - -pi(r2, ri)pi(ri, r2) 



and this function is the diagonal elements of the spinless second-order density matrix of the 
determinantal state |pi). From the second term, we have 



rWp. = ^Xj{r)xK^^){Eij}p^, Pi — > Ps 



(A6) 



and this operator is the one-body portion of the density operator, where the density operator 
is given by 



f(r) = $^X,(r)x*(r)4-. 



(A7) 



Note that we can write 



r(rU=r(r)-p,(r), 



(A8) 



indicating that f (r)p^ is determined by ps and not by pi; two different one-particle density 
matrices that yield the same density have the same f (r)p^. 
Returning to Eq. ()A4|) . from the third term we have 



f(ri,r2)p, = ^Xi(ri)x*(r2){£^ij}pi, 



(A9) 
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and this operator is the one-body portion of the one-particle density-matrix operator, given 
by 

f(ri,r2) = ^Xj(ri)x*(r2)4 = Pi(ri, rs) + f(ri, ra)^,. (AlO) 
And from the last term, we have 

f2(r2,ri)p, = ^'^Xj{ri)xK^i)Xi{r2)xl{r2){Eijki}pr, (All) 

ijkl 

and this operator is the two-body portion of the pair- function operator, Eq. ()A3|1 . 

To obtain a slight modification of Ke, we replace the determinantal state density ps, that 
appears in Eq. ()A4|) . with the exact density n, giving 



Ke' = y J dridr2r^2 P2ir2,ri)p^ + J J (iri(ir2 r^g^ n(r2)r(ri)p^ 

J f^ri(ir2rf2Vi(r2,ri)f (ri,r2)pi + J J dridr2 r{2^ f2(r2, ri)p^,(A12) 
and this operator can also be written as 

Vfe =Vee + J J dv^dr2 r^i (f (ri) - p,(ri)) , (A13) 

Replacing V^e by Vf^^ within the Hamiltonian operator, we have obtain a trial Hamiltonian: 
H,p, =H, + xj j dvidv2 r^^ Pc(r2) (f (ri) - p,(ri)) , (A14) 

where A is unity, but it can be permitted to be any constant value. Comparing this equation 
with Eq. (I2S1) yields Eq. 

APPENDIX B: ENERGY FUNCTIONAL USING INTERMEDIATE NORMAL- 
IZATION 

Using Eq. (|^. our energy functional E^, Eq. (PT|) . can be also be written as 

= - a / drg^Xr) (f (r) - p.(r)) . (Bl) 

By requiring our trial wave functions to satisfy intermediate normalization, 

(Pil^.p,) = 1, (B2) 
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we have 

E,[p^] = {p,\H,p,\^>,,,) -\j dvg,Xr) (f (r) - p,(r)) . (B3) 
This form suggest the following partitioning: 

E,[pi] = £i[piM + ^co[pi,t^] - dvg,S^) (f(r) - p,(r)) , (B4) 
where £"00 is the correlation-energy (functional) of the trial wave function: 

E^[P,M-{PAH,P.\^%;), (B5) 
and the correlation function ^^^^ is defined by 

\^.,,) = \pi) + \^%,). (B6) 
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